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Abstract. We present near-IR speckle images of 21 massive Young Stellar Objects (YSOs) associated with outflows. The aim 
of this study is to search for sub-arcsecond reflection nebulae associated with the outflow cavity. We find that 6 of the massive 
YSOs show a conical nebula which can be interpreted in terms of reflected light from the dusty walls of the outflow cavity. In 
all cases, the small scale structures seen in our images are compared with outflow indicators found in the literature. No clear 
correlation is found between the presence of the reflection nebulosity and any property such as degree of embeddedness. We 
also note that 3 of the sources show close companions, one of them belonging also to the sample with conical nebula. 
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1. Introduction 

Outflow is invariably associated with star formation and is 
likely to play a key role in the process. It can take the form 
of highly collimat ed jets that may help solve the angular mo- 



mentum problem (Bacciotti et al 



could also be present (Leeetal 
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A wide-angle wind 
and these phenom- 
ena then drive the bipolar molecular outflows in the ambi- 
ent cloud. Magnetohydrodynamic mechanisms have been pro- 
posed to explain the driving and co llimation of these out- 
flows (e.g. lGalli & Shull993llTbmisakal2002l) . This picture has 
mostly been developed for low-mass star formation. 

There are good reasons for supposing that this may not ap- 
ply to the formation of the more massive OB stars. Although 
bipolar molecular flows are just as common, there are very 
few cases where highly collimated jets are seen. On the con- 
trary, high resolution radio imaging has revealed a few cases 
where ionised winds are equatorial in n ature, i.e. perpen- 
dicular to the bipolar m olecular outflow (Ho are et al.M 1994: 
iHoare & MuxloJl99rtlHoarel2002l 

From a theoretical point of view, we may also expect 
that radiation pressure from these luminous, hot stars may 
play a larger role than magnetic pressure. Indeed, radiation 
pressure can naturally account for an equatorial win d be- 
ing d riven off the surface of an accretion disc JProga et alJ 
1998). Hydrodynamic collimation could provide an alternative 
to magnetic model s for the larger scale bipolar mol ecular flows 
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(e.g. lYorke & Welzll996UMellema & Frankll997l) . 
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One approach to achieve a better insight into these mech- 
anisms is to study the circumstellar matter distribution at the 
scales where the outflow is likely to be driven and collimated. 
In particular, the morphology of the cavity opened by the out- 
flow at scales of a few 100 AU may yield hints about the out- 
flow driving mechanism. For instance, a wider cavity would be 
expected if the outflow is rather equatorial at the base, while an 
initially jet-like outflow would produce a narrower cavity. 

A better knowledge of the circumstellar density distribu- 
tion morphology can be inferred from the near-IR observa- 
tions of the sub-arcsecond reflection nebula produced by light 
from the young star that is scattered off the dust in the cavity 
walls. These nebulae have already been observed in a few mas- 
sive YSOs f MonR2.lKoresko et al.ll993HPre bisch et alJ2002t 
S140IRS1, IHoare et alJll996t LSchertl et alJl2000t) . However 
no study has covered a large sample of sources. Observations 
of the sub-arcsecond morphology over a large sample allows 
the questions of how frequent reflection nebulae and close 
companions are amongst the massive YSOs to be addressed. 
The observed reflection nebulae can be compared with pre- 
dictions from radiative transfer models of the scattered light 
to constrain some of the properties of the circumstellar den- 
sity distributi on, e.g. optical depth cavity shape and incli- 
nation angle ( Whitnev & Hartmann 1992; Fischer & Henning 
19951 Lucas & Rochell 19981 
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tions of the ionised wind morphology from radio continuum 
observations at similarly high resolution. 

Here, we present near-IR multi-wavelength speckle ob- 
servations of the sub-arcsecond morphology on a sample of 
21 massive YSOs. The targets are known to be associated 
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with large scale bipolar outflows from molecular line and 
shocked molecular hydrogen studies. They are nearby and 
bright enough for the speckle method to b e used. Some prelim- 
inary results of this work were shown in iHoare et alJ IT996). 
In section|3the observations and data reduction are described. 
The results for individual sources are presented in Section 
In Section|4]the statistical properties of the sample (frequency 
of reflection nebulae, frequency of binaries) are discussed. 
Finally, some concluding remarks are presented in Section|5] 

2. Observations 

2.1. UKIRT Observations 

These observations (see Tabled were performed on December 
16, and 17, 1997 with IRCAM3 at UKIRT The 256 x 256 InSb 
array was used with the x5 magnifier giving a pixel scale of 
ft' 05 7 pix~'. Due to the fast readout required, for most of the 
sources, only the central 128 x 128 pixels where used, covering 
a field of view of 7'.' 3 x 7'.' 3. The seeing conditions in the near- 
IR varied between (J.' 6 and 1" 1 during the run. A minimum ex- 
posure time of 0.02 sec was used for the brighter objects, and 
up to 0.25 sec was used for the fainter ones, which is not short 
enough to completely freeze the seeing, but still gives some 
high resolution information. Observations of the target were 
alternated, normally every 500 speckle snapshots, with similar 
500 snapshots of a nearby reference star (within 1° of the sci- 
ence object) and nearby sky. More than one reference star was 
used where possible to allow for the possibility of undiscov- 
ered multiplicity in them. The log of the observations is shown 
in Table|2] 

The images were carefully cleaned of bad pixels, since they 
can contaminate the visibility function with bright peaks at spe- 
cific spatial frequencies. An average sky image (produced by 
averaging all sky frames corresponding to one target) was sub- 
tracted from every snapshot of the object and reference star. 
Cleaned and sky subtracted images were Fourier transformed 
and used to ca lculate th e powe r spectrum and the bispectrum 
(e.g. Hofmann & Weigelt 1986). The visibilities for each target 
were calculated using the corresponding reference star listed 
in Table |2] The phases were recovered with a least-squares 
method using 80 bispectral planes in the 128 x 128 image s and 
28 bi spectral planes in the 256 x 256 images jGlindemann et alJ 
Il99ll) . The visibilities were filtered using a low-pass Gaussian 
window to reduce the noise level at high spatial frequencies. 
Real images are obtained from the inverse Fourier Transform 
of the modulated visibility and the phase. The modulation of 
the visibility results in the resolution of the final images be- 
ing lower than the full diffraction limit of the telescope, but 
greatly improves the detection level of faint extended nebulos- 
ity. Generally, the reconstructed image shows artifacts, which 
can be the result of the loss of part of the speckle cloud off the 
edge of the detector, or due to a poor phase recovery. They usu- 
ally appear concentrated close to a strong peak, and are at the 
level of a few percent of the brightest peak. Tests on reference 
stars were used to check for the level of these artifacts, as well 
as for binarity amongst the reference stars themselves. 



2.2. CalarAlto Observations 

The observations were made at the 3.5m telescope on Calar 
Alto, on December 9, 10 and 1 1, 1994. The MAGIC 256x256 
NICMOS 3 camera was used at a pixel scale of 0'.'073 per pixel 
in order to sample the diffraction limit at K. For simple point 
sources 1500 128x256 frames with an integration time short 
enough to freeze the seeing were accumulated on the object in- 
terspersed with a similar number on one or more reference stars 
less than 1° away. The seeing was approximately 1 arcsec. Only 
half of the array was saved in order to save disk space and speed 
up the readout. The source was alternately placed in lower and 
upper quadrants, the blank one being used for sky subtraction. 
Where multiple sources or complex structure extended across 
the field of view the whole array was saved and separate sky 
frames where obtained. The reduction procedure is similar to 
the one used for the UKIRT images. 

A typical value for the limiting surface brightness in our 
observations can be obtained from the seeing-limited magni- 
tudes listed in Tableland the typical values for the resolution 
and noise level in the restored images (see figures fTlto ll5b . As 
an illustration, we take the instance of GL 490 in the H band 
from our Calar Alto campaign, since the reconstructed image 
(see Fig. |8j is dominated by a strong point source with only 
a faint nebulosity and no other sources present in the field of 
view. With a seeing-limited magnitude of 7.6 and a restored 
resolution of ff.'23, only 8% of the flux is in the central pixel. A 
lcr noise level of 2.5% of the brightness peak is found within 
1" of the peak, whilst %\% noise is found further out. These 
yield 3<x limiting surface brightnesses of 7.4 mag/arcsec 2 and 
10.1 mag/arcsec 2 respectively. These estimates are similar for 
the rest of the sources for which we could reconstruct a high 
resolution image. Therefore, the limiting surface brightness (in 
units of mag/arcsec 2 ) in the region of interest (i.e. within 1" of 
the central peak) is of the same order as the total seeing mag- 
nitude (in units of mag). This is less valid for sources with a 
strong nebulosity or for fields with multiple sources. The noise 
level is determined by the dynamic range (~ 100) of the speckle 
technique used here. 

3. Results 

In this section we present the results on individual objects be- 
fore discussing the sample as a whole in the next section. 

3.1. MonR2IRS3 

Mon R2 IRS3 is the brightest of a cluster of compact IR sources 
with associated extended emission located at a distance of 
950 pc. IRS3 is located approximately 3 0" to the east of a 
blister type Hn region of ~ 27" diameter (Ma ssi et"afl [T985). 
The cluster is located near the centre of a giant CO outflow 
(Ball y & Ladall983UWolf et alll9 90) oriented at a PA ~ 135° 
with the approaching gas flowing to the NW. Due to the low 
resolution of the CO observations, it is not clear which of the 
IR sources is powering this outflow. On scales of a few arcsec- 
onds, the IR reflection nebula is extended to the south and to 
the east. The polarisation pattern shows centro-symmetry in- 



Table 1. Target list. 
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UxSJliL. 1 


KA 


UrA 


L(L G ) 


Kei 


D (KpC ) 


Kei 


mj 


™ d 

m H 


m K 


GL4029 IRS1 


03 01 31.3 


+60 29 13 


2.1 x 


10 4 


15,17 


2.20 


16 


12.9 


10.0 


7.7 


GL437 


03 07 24.6 


+58 30 44 


3.0 x 


10 4 


18 


2.70 


5 


12.9 


11.3 


9.8 


GL490 


03 27 38.8 


+58 46 60 


1.4 x 


10 s 


7 


0.90 


19 


10.0 


7.6 


5.3 


LkHor 101 


04 30 14.4 


+35 16 25 


1.2 x 


10 4 


7 


0.80 


6 


8.1 


5.7 


3.1 


S235IRS4 


05 40 52.3 


+35 41 31 


1.2 x 


10 3 


10 


1.80 


9 


11.0 


9.9 


9.0 


S235IRS1 


05 41 06.9 


+35 49 39 


1.5 x 


10 3 


10,11 


1.80 


9 


12.2 


10.2 


8.1 


S235 IRS2 


05 41 11.0 


+35 50 02 


2.5 x 


10 3 


10 


1.80 


9 


12.2 


8.8 


6.5 


NGC 2024 IRS2 


05 41 45.8 


-01 54 29 


1.1 X 


10 4 


21 


0.40 


20 


11.6 


7.4 


4.6 


MonR2IRS2 


06 07 45.8 


-06 22 54 


5.0 x 


to 3 


22 


0.95 


1,2 


18.0 


14.1 


9.5 


MonR2IRS3 


06 07 48.1 


-06 22 55 


6.0 x 


10 3 


1 


0.95 


1,2 


12.7 


9.5 


6.6 


GL5180 IRS1 


06 08 53.4 


+21 38 29 


1.1 X 


10 4 


23 


1.50 


23 


11.6 


10.7 


10.3 


S255IRS1 


06 12 53.9 


+ 17 59 23 


3.2 x 


10 4 


25 


2.40 


24 


16.5 


12.9 


9.6 


S255 IRS3 


06 12 54.0 


+ 17 59 23 


1 9 v 
J.Z X 


i n 4 
1 u 


25 


2.40 


24 


11.3 


10.6 


10.4 


GL961-W 


06 34 37.4 


+04 12 43 


6.3 x 


10 3 


26 


1.60 


8 


10.4 


8.9 


7.7 


GL961-E 


06 34 37.7 


+04 12 44 


1.5 x 


10 3 


26 


1.60 


8 


14.6 


10.5 


7.8 


GL 961 -Fan 


06 34 35.6 


+04 12 46 








1.60 


8 


13.1 


11.4 


9.7 


GL989 


06 41 09.7 


+09 29 36 


4.0 x 


10 3 


27 


0.75 


14 


11.0 


7.6 


4.8 


S106IR 


20 27 26.8 


+37 22 48 


2.0 x 


10 4 


29 


0.60 


28 


10.4 


7.7 


5.9 


GL2591 


20 29 24.9 


+40 1 1 20 


4.0 x 


10 4 


30 


1.50 


30 


14.4 


11.0 


6.1 


S 140 IRS 1 


22 19 18.1 


+63 18 47 


5.0 x 


10 3 


3 


0.90 


4 


11.7 


8.6 


6.2 


NGC 7538 IRS 1 


23 13 45.4 


+61 28 12 


5.0 x 


10 4 


12 


2.70 


13 


15.3 


14.1 


8.9 



a Equinox J2000. 

b References for the luminosity. 

c References for the distance. 

d Magnitudes extracted from the literature. 



Cohen & KuhiN 19771) : 6. 


HerbisUl97l|) : 7 lHarvev et alj<1979b|) : 8.lBlitz & Thaddeusl 


Jl98d) : 9.lEvans & BlairlJl98ll) : lO.lEvans et al 


jl981 


) ; 11. Nordh et al. 


1 1984) : 12. IcamDbell & ThoniDsorJ Il984l) : 13. iBlitz et al. 


1 1983) : 14. IPark etalJ feOOCl) : 15. ISnell et al 


(1988 


) ; 16. Becker & Fenkart 1 1971) ; 17. Thronson et al. 1 1980b) ; 18. Wvnn- Williams et al. (1981) ; 19. Harvev et al. 1 1979a) ; 20. 


Anthonv-Twarojl 1982) : 


21.|jiang et alJ<1984l) : 22 iHoward et . alHl994l) : 23.IWu et al 


J<1996h : 24.lltoh et alJfcOOll) : 25.lHoward et all 



( 1997) ; 26. Castelaz et al. 1 1985) ; 27. Harvev et al. ( 1977) ; 28. Staude et al. 1 1982) ; 29. Harvev et al. 1 1982) ; 30. Wvnn- Williams 1 1982) 



dicating that IRS3 is the illuminating source. It als o shows a 
polar isation disc elongated in the SE-NW direction dYao et alJ 
119971) . 

Thi s object has been inve stigated with speckle imaging be- 
fore by Kores ko et al.l (11993) at 2.2, 3. 8 and 4.8 /urn and more 
recently in deep images at H and K by Preibisch et al. (2002) 
using the 6 m SAO telescope. We present our K and H-band 
speckle images in figures ^ ar, d ID f° r comparison with the 
rest of our survey. In the K-band two stars are clearly present, 
IRS3 S and IRS3 N, that are s eparated by 0787 + 0703, in good 
agreem ent with the works o f iMcCarthvl ill 9821) . iKoresko et alJ 
i ll 993b and IPreibisch etatl J2002I) . IR S3N is actually com - 
posed of two stars (sources B and C in Preibis ch et alJ l2002). 
These two stars are also seen by Koresko et al. in the K band. 
However, they are only marginally resolved in our K band 
data. The conic al reflection nebula e manating to the S-SE of 
IRS3 S found by Kores koet alJ 1 19931) is clearly visible in both 
the H and K bands. At H, the brightest of the two stars in 
IRS3N (source B in the work by Preibisch et al.) disappears 
completely, whilst the star illuminating the southern nebula is 
faint, if present at all, at the northern tip of the conical neb- 
ula. These monopolar nebulae are interpreted as light scattered 



off the walls of an outflow cavity that is directed somewhat 
towards us. The other lobe that would be associated with the 
redshifted outflow is not visible due to high extinction in an in- 
ferred flattened envelope perpendicular to the bipolar outflow 
direction. Figure 13.11 represents a colour-coded image, which 
stress the presence of the highly reddened stars and the rela- 
tively blue conical nebula associated with IRS3 S. The scatter- 
ing efficiency of dust increases rapidly with decreasing wave- 
length through the near-IR, causing the reflection nebula to be 
relatively blue. 

The morphology of the southern nebula indicates that 
in this case the outflow is at a position angle (PA) of 
198° + 15° wit h the a pproaching gas towards the south-west. 
IPreibisch et ail d2002l) show that the emission to the E of 
IRS3 N is the combination of another star and a jet-like feature 
emerging from IRS3 N at a PA ~ 50°. Neither of these outflow 
orientations correlate w i th the large scale molecular outflow 
described bv lWolf etall dl990l) (PA ~ 135° with the blue gas 
moving towards the north-west). Therefore, the IRS3 sources 
are not likely to be driving the large scale bipolar molecular 
outflow in the region and we need to look for other outflow in- 
dicators nearer to IRS 3 itself. One should keep in mind though, 
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r liter / 


Standard 


Average 




Number 


Notes 




Telescope' 1 


otai 


counts 


Frames 




GL4029 IRS1 


K/C 


SAO 12585 /SAO 23728 


120 


0.20 


1536 


N! 


GL437-S 


K/U 


CMC 302423 / CMC 601567 


21 


0.13 


1200 


N 


GL490 


I/U 


CMC 401688 / CMC 302688 


31 


0.20 


1500 


_ 




H/U 


SAO 24017 /SAO 23967 


116 


0.05 


4000 


N! 




K/U 


SAO 24017 


474 


0.02 


3500 


U 




Bry/U 


SAO 24017 


107 


0.13 


2000 


U 




Con/U 


SAO 24017 


70 


0.13 


2000 


U 




H/C 


SAO 23999 c / SAO 24017 


120 


0.06 


1280 


N 




K/C 


SAO 24098 / SAO 24017 


537 


0.06 


1000 


U 


LkHff 101 


I/U 


SAO 57254 / SAO 57264 


73 


0.05 


3000 


C 




H/C 


SAO 57239 / SAO 57264 


692 


0.06 


1536 


c 


S235 IRS4 


H/U 


SAO 58357 / SAO 58383 


65 


0.13 


3000 


u 




K/C 


SAO 58343 / SAO 58272 


54 


0.06 


1536 


u 


S235IRS1 


H/UC 


SAO 58357 


56 


0.20 


2000 


N 




K/U 


SAO 58343 


119 


0.05 


2000 


N 




K/C 


SAO 58343 / SAO 58272 


122 


0.20 


1536 


N 


S235 IRS2 


K/C 


SAO 58343 / SAO 58272 


190 


0.06 


1536 


U 




H/C 


SAO 58357 


107 


0.10 


1536 


u 


NGC 2024 IRS2 


K d /C 


SAO 132480 


853 


0.06 


1536 


u 




H/C 


SAO 132410 / SAO 132434 


144 


0.06 


1536 


u 


MonR2IRS2 


K/U 


SAO 132881 / SAO 132874 


15 


0.20 


1536 


u 


MonR2IRS3 


K/C 


SAO 133884 / SAO 132874 f 


106 


0.06 


1536 


N 




H/U 


SAO 132874 


139 


0.20 


2000 


N 




K/U 


SAO 132884 


75 


0.02 


2000 


N 


GL5180IRS1 


K/C 


SAO 78056° / SAO 77997 


23 


0.3 


768 


U 


S255IRS1/3 


K/C 


SAO 95350 


41 


0.20 


1536 


N 


GL961-E/W 


H/U 


SAO 1 14047 /SAO 11 4094 


209 


0.25 


1000 


N!/C 




K/C 


SAO 1 14047 / SAO 1 14O80 


62 


0.06 


1792 


U 


GL 961 -Fan 


K/U 


SAO 1 14047 / SAO 1 14094 


88 


0.13 


2000 


u 


GL989 


H/U 


SAO 95987 


140 


0.05 


3000 


u 




H/C 


SAO 95987 


120 


0.06 


2048 


N! 




K/C 


SAO 114234 /SAO 114245 


68 


0.06 


1792 


U 


S106IR 


K/C 


SAO 76077 / SAO 70120 


330 


0.06 


2048 


U 




H/C 


SAO 70 120 /SAO 70 129 


109 


0.13 


1280 


u 


GL2591 


H/U 


SAO 70099 / SAO 49777 


29 


0.13 


2000 






K/U 


SAO 70052 


381 


0.05 


1500 


N! 




K/C 


SAO 70052 / SAO 49730 


274 


0.06 


2048 


N! 


S 140 IRS 1 


K/C 


SAO 20051 /SAO 19948 


250 


0.06 


1600 


N 


NGC 7538 IRS 1 


K/C 


SAO 20455 / SAO 20563 


21 


0.20 


1536 


N! 



a U: UKIRT. C: 3.5 m telescope at Calar Alto. 
b Average peak counts per frame. 

c These stars were found to be close binaries and so not included in the reduction. 
d Object was too bright at K and so a narrow-band filter at 2.26yum was used instead. 

c Main features observed in the reconstructed images. N: Bright nebula likely to be associated to the outflow cavity. N!: Very low level or 
hardly resolved nebulosity. C: Close companion (with a separation < 1'.'5). U: Unresolved source. -: Reconstruction not possible due to the 
low number of photons per individual frame. 



that single-dish CO observations have a very low spatial reso- 
lution (several arcseconds) compared with our speckle images. 
The giant molecular outflow in the Mon R2 region maybe a su- 
perposition of several smaller outflows with slightly different 
orientations, each one driven by a different source. Therefore, 
any attempt to relate the near-IR sub-arcsecond morphology 
with large scale structures should be taken with caution, even 



though this is an interesting exercise to search for possible cor- 
relations in orientation. 

It is also interesting to compare the orientation of the 
cavit y with the distribution of methanol masers presented 
by iMinier etal . (2000). They show a group of masers form- 
ing a linear structure with an approximate length of (X'17 at 
PA = 230° + 20°. They interpret this structure in terms of a 
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Mon R2 IRS 3 (K band) 

1" 




N A 

-4 

E 



Fig. 1. Reconstructed image of Mon R2 IRS3 in the K band. A 
Gaussian window of FWHM=5.3 arcsec -1 was applied to the 
visibility to achieve a resolution of 07 19. The contours are at 5, 
10, 15, 30 and 50 % of the peak. The noise level is ~ 0.4 %. 
The sources are labelled using the notation by Koresko et al. 
1993. The arrows indicate faint stars used to align the H and 
K images. In this and in the following reconstructed images 
the intensity scale is normalised to the brightest peak with the 
brighter regions represented by the darker greys. 



Mon R2 IRS 3 (H band) 

& 




N A 

-4 

E 



Fig. 2. Reconstructed image of MonR2IRS3 in the H band. 
The resolution is 0719. The contours are at 5, 10, 15, 30 and 
50 % of the peak and the noise level is ~ 1 %. The labels used 
are the same as for the K band image. The star symbols indi- 
cate the position of the stars IRS3 N and IRS3 S in the K band 
image. The small arrows indicate the position of the two stars 
used to align this image and the K band image. 



nearly edge-on circumstellar disc. This would imply a sys- 
tem with a rotation axis in the NW-SE direction, which is in 
complete disagreement with what the near-IR sub-arcsecond 
nebula suggests. Alternativ ely, the situation m aybe similar to 
other sources (e.g. W75N Torrel les et alJ ll997) where the lin- 
ear maser morphology appears to trace the outflow axis. 

Final ly, the ori entation of the polarisation structure shown 
bv lYao et alJ Jl997t) is considered. Their seeing-limited images 
indicate the presence of a polarisation disc at PA ~ 140°, which 
is roughly perpendicular to the maser structure. It is likely that 
IRS3 S, which appears to be the strongest source of scattered 
light in our UKIRT speckle images, is the main contributor to 
the polarisation disc. However, high resolution polarimetry in 
the near-IR is necessary to resolve the polarisation structure of 
each of the components independently. 



Mon R2 IRS 3 (H + K) 




3.2. S140IRS1 

S 140 IRS 1 is located near the centre of a bipolar molecular out- 
flow at a PA ~ 160° with the blue- shifted lobe pointin g towards 



the SW JBlair et al.1 


19781 


Ballv & Ladalll983l lHavashi et alJ 


ll987llMinchinetalJ 


1993 




Minchin et al. 


(1995 


find that the 



C O emission at the ambient cloud velocity traces the cavity 
wall associated with the blue-shifted CO outflow lobe. 

Seeing-limited imaging polarimetry in the optical and in 
the IR show that S 140 IRS 1 is the brightest member of a clus- 



Fig. 3. Colour-coded image of MonR2IRS3. The emission at 
H and K are represented by the blue and red colours respec- 
tively. The blue colour of the nebulae indicates the scattered 
nature of the light. 
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ter of embedd ed sources surrounded by an exten ded reflection 
nebula jForrest & Shurejl98rllLenzerJl987tlEvans et alJ l989: 
lHarker et al.lll997HYao et al.lll998h . The polarisation structure 
of the nebula indicates that IRS1 and IRS3 are the main illu- 
minating sources of the nebula, but the polarisation structure is 
complex. 

Radio continuum observations reveal a thermal source 
with a spectral index of 0.8 indicative of an ionised wind. 
Observations that resolve the source show that the radio emis- 
sion is elongated perpendicular to the larg e scale outflow di- 
rection ( lSchwarta ll989: Tofani et al. 1995). The higher reso- 
lution observations bv lHoare & Muxlowl i 1996) show a struc- 
ture at a PA of 44° ±1°. Multi-epoch observations show that 
this st ructure is not a radio jet as first propose d bv ISchwartzl 
Jl989t) . but an equ atorial wind jHoard as was also 

found for S106IR faoare et al.lll994l) and possibly GL490 
JCampbelletalJl986l) . ~ 

Our K-band speckl e image of S 140 IRS 1 is shown in fig- 
ure|4](see also lHoare et alJ 19961) . This image shows a monopo- 
lar nebula that emerges from IRS1 and extends towards the 
SE at a PA ~ 150°+ 15°. The azimuthally averaged visibility 
function indicates that 22 ± 2% of the flux arises from the 
point source, whic h is likely to be the central proto-star (see 
ISchertl etalj EoOO). The morphology in figure [4] agrees with 
deeper speckle imaging bv lSchertl et alJ (|2000). Their polari- 
sation data also show the centro-symmetric pattern that con- 
firms that this is a reflection nebula illuminated by S 140 IRS 1. 
This is consistent with the light being reflected off the walls of 
the blue-shifted outflow cavity in a classic monopolar reflec- 
tion nebula. All this fits in with the interpretation of the radio 
structure as an equatorial wind. 

The only questio n mark over this picture is the discovery by 
Wei geltetalJ(l2002h of a set of three bow shock-like structures 
to the NE of S 140 IRS 1 in wide-field K-band speckle images. 
They interpret these as arising from precessing, episodic jet ac- 
tivity from S 140 IRS 1. This would require a return to the jet 
interpretation for the radio emission, which is inconsistent with 
the direction of the main outflow indicators, i.e. the CO outflow 
and reflection nebula. One solution to this conflict would be if 
the bow shock structures were driven by a jet from a lower- 
mass YSO located near to S 140 IRS 1 . 



3.3. S255IRS1 and IRS3 

S255 IRS 1 and IRS3 form part of a cluster of IR sources located 
between the evolved Hn regions S255 and S257 at a distance of 
2.5 kpc. A molecular outflow along the N-S direction with its 
approaching l obe towards the sou th appears to be centred near 
the IR cluster jHeveretalJl9 89). Infrared polarisation images 
bv lTamuraet alJ dl99lh show a bipolar nebula (IRN1) at PA of 
63° illuminated by IRS3, which indicates the presence of an 
outflow in the this direction. They also find another reflection 
nebula (IRN2) to the north of IRS 1 that is illuminat ed by that 
source. Thermal IR observations bv lltoh et alJj200l|) show that 
IRS3 is the more luminous and embedded object. Miralles et al. 
il 19971) find molecular hydrogen emission from HH-like objects 
along the jet direction. 




Fig. 4. Speckle reconstructed image of S 140 IRS 1 in the K band 
at a resolution of 0'.'19. The contour levels are at 1%, 2.5%, 4%, 
5.5%, 7%, 8.5% and 50% of the peak. The noise level is < 1%. 
The grey-scale varies from -2% (white) to 8% (black). 



The K-band reconstructed image of this region is shown in 
figure |3 The field of view contains the point sources IRS1, 
IRS3, IRS22 and the two refl ection nebulosities I RN1 and 
IRN2 using the nomenclature of Tam ura et alJ lll99lh . IRN1 is 
seen to be very conical in shape with an opening angle of 35° 
and at PA~ 58° from its illuminating source IRS3. Although 
the main orientation of the CO outfl ow is along the N-S di- 
rection, the CO observations of lHever et alJ dl989h show some 
evidence for high velocity blue-shifted gas towards the SW 
roughly along the outflow axis indicated by IRN1 in figure|5] 

IRN2, which Tamura et al. deduced was present from the 
polarisation pattern, is clearly resolved in our observations. 
There is a gap between the main IRN2 nebula and IRS1 and 
its morphology is not entirely consistent with IRS 1 being its il- 
luminating source. The presence of foreground clumps or intra- 
cluster matter could be a possible explanation for this gap. 
Significant artifacts exist on the RA axis through IRS 1 so that 
the reality of features close to IRS1 are questionable. Some 
spatially distinct diffuse emission is seen on the opposite side 
of IRS 1, which maybe the other reflection lobe. If these nebu- 
losities were associated with the main N-S outflow, we would 
expect the southern blueshifted lobe to be brighter. 

3.4. GL437-S 

GL437 is a compac t cluster of IR sources 
JWvnn-Williams et all 19 8 1 ) located a distance of 2.7 kpc. The 
cluster, with a total luminosity of ~ 3 x 10 4 L , comprises 
four bright sources denominated GL437-N, S, E and W. The 
cluster is located at the centre of a poorly collimated molecular 
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Fig. 5. Speckle reconstructed image of S255 IRS 1 in the K band Fig. 6. Speckle reconstructed image of GL 437-S in the K band. 



at a resolution of 0729. The contour levels are at 0.5%, 1%, 
1.5% and 2% of the peak. The noise level is < 1%. The grey- 
scale varies from -1% (white) to 4% (black). The negative fea- 
ture (-3.5%) west of IRS 1 is a 'ghost' of IRS3. 



The resolution attained is 0722. The colour scale varies from - 
1% (white) to 5% (black) of the brightness peak. The contours 
are at 15%, 20%, 25%, 30% and 50% of the peak, and the noise 
is at the 4% level. 



outflow JGomez et all 1 9921) . Iweintraub & KastneJ i ll 9961) find 
that the cluster is embedded in an infrared polarised nebula. 
The polarisation pattern is mainly centro-symmetric with 
respect to the embedded source WK 34, which is close to 
GL437-N. How ever, centro-symmetry can also be observed 
about GL 437-S dWeintraub et alll996l) . 

Figure |6] shows the reconstructed image of GL 437-S at a 
resolution of 0722. The data indicate the presence of a sub- 
arcsecond monopolar nebula extended towards the SE at a 
PA = 135° + 10°. The noise level is high in this image since 
GL 437-S is quite faint, but the extended emission is bright 
relative to the star. However, the presence of the nebulosity is 
consisten t with the centro-sym metry in the near-IR polarisation 
seen by ( Weintrau b et al.lll996l) . The high-velocity CO emis- 
sion shows several peaks (see figure 3 of Gomez et al. 1992). 
Although the main outflow axis appears to be in the N-S di- 
rection, the presence of multiple outflows in the region cannot 
be discarded. In particular, there is some blue-shifted emission 
towards the SE, and red-shifted gas to the NW, that may trace 
a bipolar outflow oriented roughly along the same direction as 
the sub-arcsecond nebula shown in figure|6] 

3.5. LkHa 101 



LkHa 101 is a Herbig Be star with strong Ha emission (Herbig 
1 19561) which illuminates the reflection nebula NGC 1579 
jRedmanetalJll986l). It is the brightest of a cluster of IR 
sources ( Bars onv et al.ll99ll) at a distance of 800 pc. Its spec- 
tral type corresponds to a late O or early B star JBrown et alJ 



ll976l | Harvev et alJl979bll. It excites a weak Hh region of size 
~ 1' Becker & Whitelll988l) . LkHa 101 is also the source of 
a compact ionised wind. The radio continuum emission from 
the wind has been studied at a wide range of frequencies 
dBrown et alJl976HBieging et alJl984lBecker & Whitell 988) 



showing a 



„0.6 



law. Hoare & Garrington ( 1995) resolve the 



wind into a clumpy and approximately circular distribution on 
the sky. 



Ballv & Lada ( 1983) studied the large scale molecular gas 
motions in the region. They classified LkHa 101 as an inter- 
mediate velocity source (with velocities in the 10-30 kms _I 
range). However, the existence of a bipolar molecular outflow 
was not cle ar from their data, and has been completely dis- 
carded by iBarsonv et alJ Jl990l) . This object shares many of 
the common features of massive young stellar objects, but the 
fact that it is seen in the optical indicates that it is probably in 
a fairly advanced evolutionary stage. 

The sub-arcsecond structure in the near-IR of LkHa 101 
has been studie d by several authors using ID speckle inter - 
ferometry dPewarf & Dvcklll993l Leinert et alJll997l l200ll) . 
They classify the source as either unresolved or with a core- 
halo morphology. However, this object was resolved into two 
components (the brightest surround ed by a disk) with a sepa- 
ration of 07 1 8 and at a PA ~ 70° by Tuthil l et all d200ll 120021) 
using aperture-mask interferometry with the Keck I telescope 
at a resolution of ~ 07023 . 

For completeness, we show in figure 7 our J band speckle 
data of LkHa 101 taken at UKIRT at a resolution of 07 1 1 . A 
separation of 0717 ± 0702 between both components was mea- 
sured directly in the reconstructed image and also by fitting the 
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Fig. 7. Reconstructed image of LkHalOl in the J band 
with a resolution of 0711. A Gaussian window of 
FWHM=9.1 arcsec -1 was applied to the visibility to re- 
duce the noise associated with the high spatial frequencies. 
The faint component is located 0/17 away (~ 130 AU at a 
distance of 800 pc) and it is ~ 5 times fainter than the primary. 
The background noise is ^ 1% of the peak. 



bispectral phase. The position angle of the fainter component 
is 75° + 5°. Aperture photometry on the reconstructed image as 
well as fitting the bispectrum yielded a brightness of the fainter 
component with respect to the total of 0.17 + 0.03. The binary 
system can also be seen in the H band data taken during the 
Calar Alto campaign. A separation of 0/19 + O/03 and a rela- 
tive brightness of 0.06 + 0.02 with respect to the total. These 
values are consistent with those measured by Tuthill et al. at 
higher resolution. 

3.6. GL490 

This massive YSO is associated with a parsec scale mo lecular 
outflow that was first detected bv lLada & Harve vl(ll98ll). Mor e 
detailed structure of the out flow is given bv lSnell et al.l Jl984). 
iMitchell et all 119921 17995). All observations agree in that the 
outflow axis lies in the NE-SW direction (PA ~ 45°) with the 
blue-shifted lobe towards the south-west. The CO outflow is 
not highly collimated, with a considerable overlap of the blue- 
shifted and red-shifted gas. 

An stellar-like optical counterpart is located at the apex 
of a conical reflection nebula in the blueshifted outflow lobe 
dCampbell et alJfT986l) . The light centre of the optical source 
was found to shift to the SW at bluer wavelengths consistent 
with a scattered rathe r than direct source for the optical light 
(Campbe ll et alll98 6). Near-IR imaging polarimetry of the re- 



gion performed bv lMinchin et alJ dl99lh shows a cometary re- 
flection nebula extending for about 20" towards the south-west, 
with its head at the location of GL490. This nebula shows a 
wider opening angle than its optical counterpart, and it is char- 
acterised by a centro-symmetric polarisation pattern with a po- 
larisation disc at a PA of 120° ± 2°. 

An ionised stellar w ind is seen in the radio with a 
spect r al index of ~ 1 .2 dSimon et alJ Il983l ICampbell et alJ 
19861 iHenninget all [l990). The radio source appears ex- 
tended i n a direction perpe ndicular to the outflow at 
15 GHz JCampbell etalj fl986). At the position of the star, 
iKawabe et afl ]l984) found CS emission elongated in a direc- 
tion roughly perpendicular to the outflow axis that can be inter- 
preted as produ ced in a dense d isc or torus of material around 
the IR source. Schrever et al ] 120021) also find an elongated 
structure at a PA~-45° in the CS (J=2-l) transition. Further 
evidence of the disc is obta ined from the 2.7 m m continuum 
emission detected by M undv & Adelman1 ( ll988l) . which is also 
elongated perpendicular to the outflow. 

The sub-ar csecond structure of GL 490 in the near-IR was 
first studied by Ha as et alJdl992i> . Their ID speckle interferom- 
etry in the H, K and L bands shows that the source is formed 
by an unresolved core of ~ 07 1 and an extended halo of ~ 172. 
The halo appears relatively bluer, slightly elongated towards 
SW (PA~ 125°), and shows a centro-symmetric polarisation 
pattern, which indicates a scattered ori gin for the emission. T he 
same core-halo structure is found bv lDewarf & Dvckl dl993h at 
3.8 /im. 

Our K-band speckle imaging does not resolve any sub- 
arcsecond structure. However, extended emission is seen in 
our H band reconstructed image taken at Calar Alto (see fig- 
ure |8j. This blue colour of the nebula strongly indicates that 
it is scattered light. The H-band nebulosity is mainly perpen- 
dicular to the outflow (PA~ 120°), but appears to curve off 
towards the blueshifted SW lobe at the ends. Such a morphol- 
ogy is consistent with arising at the base of the wide opening 
angle blueshifted outflow cavity. This structure is very similar 
in size and orientation to the extended radio emission seen by 
ICampbelletaillll986l) . 

For the speckle data taken at Calar Alto (see Table |2ji, it 
turned out that one of the reference stars used (SAO 23999) 
was actually a binary and so could not be used in the recon- 
struction. This resulted in somewhat higher than normal arti- 
facts and partly for this reason the speckle observations of GL 
490 were repeated at UKIRT. Figure [9] shows the H band im- 
age reconstructed at a resolution of 0/18 from this dataset. The 
elongated morphology seen in the image taken at Calar Alto 
appears to have vanished in the new UKIRT image. Tests done 
with cross-calibrated PSF stars favour the reality of the disap- 
pearance of the nebula in GL 490. 

Interestingly, the extended radio emission observed by 
ICampbell et alJ dl986h has also not been detected in more re- 
cent deep 8 GHz VLA observations by Hoare & Gibb (in 
prep.). Hence, it is possible that these changes in the morphol- 
ogy of near-IR and radio observations maybe caused by real 
variability in the source. One possible explanation could be ob- 
scuration near the star (above the disc) preventing the sides of 
the cavity being illuminated with near-IR light, thereby reduc- 
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Fig. 8. H band speckle reconstructed image of GL490 from 
data taken at Calar Alto in 1994 The resolution is (X'23. The 
grey-scale varies between -2.5% (white) and 9.4% (black) of 
the peak. The contours are at 1.5, 3, 4.0 and 50% of the peak. 
The noise level is 0.2% away from the source. 



GL 490 (H band) 



1" 



Fig. 9. H band speckle image of GL490 taken at UKIRT 
in 1997. The resolution achieved was 0'.'18. The visibility 
was multiplied by a Gaussian window with a FWHM of 
5.5 arcsec -1 . The grey-scale varies from -1% to 3%. The con- 
tours are at 1.5, 3, 4.5 and 30% of the peak. The most negative 
feature is at a -1.5% level. The noise is ~ 0.2% away from the 
source. 



ing the reflection nebula and EUV radiation, thereby stopping 
photoionisation and hence the radio emission as well. Further 
multi-epoch observations at different wavelengths are needed 
to probe this possible variability. 

Any reflection nebulosity should be brighter relative to the 
star at J compared to H. Unfortunately, GL 490 is too faint at 
J for the reconstructed image taken at UKIRT to confirm the 
reality of the nebulosity. Conversely, the source is so bright at 
K that speckle imaging could be performed through a Bry filter. 
However, no extended emission due to ionised gas was found 
here either. 

3.7. GL961 E, W and Fan 



GL961 is a luminous (Castela z et al.lfl985l) young stellar ob- 
ject that is located in the south-west rim of the Rosette nebula 
at a distance of 1.6 kpc. This source is very faint in the optical 
feiroalll98lllBeetz et alJll976l) . and lies approxima tely 30" E 
of an optical fan-shaped nebula (GL961-FanV iLenzen et alJ 
found that GL961 is actually composed of two close 
IR objects, separated by about 5" (8000 AU) at a PA ~ 70°. 
ICastelaz et alJ i 19851) show that the spectral energy distribu- 
tion of GL 961-E dominates at wavelengths longer than 2.2 yum, 
while the western object is brighter at shorter wavelengths. 

The double source is located at the centre of a parsec 
scale molecular outflow oriented in the N-S direction with th e 
approaching lobe towards the north llLada&Gautierl U982). 
Observ ations of the H2 emission in the region performed by 
lAsninl |l998) show multiple bow-shocks emerging from the 
double system. One system is oriented at PA ~ 40° and the 
other one at PA ~ 0°, which corresponds with the orientation 
of the CO outflow. 

Figure [K)] shows the speckle reconstructed image of the 
double system formed by GL 961-E and GL961-W. In this 
case, the field of view is 14'.' 6 x 14'.'6 since the full 256 x 
256 pixel area of IRCAM3 was used. Both sources can 
be clearly identified at a separation of 5'.'29 + $.'03 and a 
PA = 17.9° + 0.3°. A very faint third source (GL961-Wb) can 
also be seen at a separation of 1'.'48 + a.' 06 and a PA = 59° + 3° 
from GL961-W (renamed GL961-Wa here). GL961-Wa con- 
tributes a 14.2% ± 0.9% to the total flux, while the contribution 
of GL961-Wb is 1.0% + 0.2%. We use the magnitudes given 
in Castel azet al.1 dl985l) to estimate a total magnitude in the H 
band of 8.65+0.01 for the system GL 961-E, -Wa and -Wb. This 
yields a H magnitude of 10.8+0.1 for GL961-Wa and 13.6+0.1 
forGL961-Wb. 

Some extended emission at the 1% level can be seen around 
GL 96 1 - Wa. The shape of the extended emission is fairly uncer- 
tain since it is very close to the noise level. In summary, there 
appears not to be any clear feature in this source that can be 
related to any of the outflow activity indicators seen at a larger 
scale. 

High resolution observations of GL 961 -Fan were also per- 
formed at K, since it may constitute an independent site of star 
formation. However, no sub-arcsecond extended emission was 
found at a level S 1 % of the brightness peak. 
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Fig. 10. Reconstructed image of GL 961 at H. A Gaussian win- 
dow with a FWHM=3.1 arcsec -1 was applied to the visibility 
function in order to filter out the noise associated with the high 
spatial frequencies. The image has a resolution of 0732. The 
contour levels are at 1%, 3%, 5%, 10%, 30% and 50% of the 
brightness peak. The background noise is ~ 0.1%. 

3.8. S235 IRS1 , 2 and 4 

These luminous YSOs dAllenlll972h are located in the S235 
molecular cloud at a distance of 1.8 kpc. S235IRS1 was ob- 
served in speckle mode during 1994 at the Calar Alto observa- 
tory in the K band, and in 1997 at UKIRT in the H and K band. 
The K band speckle reconstructed images from these two ob- 
servations at a resolution of 0727 are shown in figuresfTTh and 
II lb . In both K images, S235 IRS 1 can be seen near the cen- 
tre, and two other stars are detected to the west of IRS 1 . One 
of them is located 4'.' 13 + 0707 to the NW of S235IRS1, and 
the other 370 ± 071 to the SW. The brighter of these stars con- 
tributes a 2.2% ± 0.1% to the total flux, while the contribution 
of the fainter is 0.6% + 0.1% . Assuming a total magnitude of 
10.25±0.07 (11" aperture, c.f. lEvans & Blairll98ll) . the magni- 
tudes of the bright and faint stars are 12.29+0.08 and 13.7+0.2 
respectively. 

Some extended emission elongated in the E-W direction is 
clearly seen in S235 IRS 1 . The extended emission that is de- 
tected at 173 ± 071 towards the east of the central peak, is cer- 
tainly real. The contribution to the total flux of this feature, 
measured on the UKIRT image, is 6.2% ± 0.4%. However, the 
features closer in are more questionable since there are negative 
artifacts at a similar level. Intrinsic time variability of the neb- 
ulosity cannot be discarded, since the extended emission de- 
tected to the west of IRS1 at Calar Alto (see figure [TTV) seems 
to have vanished in the 1997 data. 

The reconstruction of the H band data yielded a low SNR 
image due to the low level of counts on the reference star. 



Nevertheless, inspection of the seeing-limited image indicates 
that the nebulosity to the east of IRS 1 is relatively bluer than 
the brightness peak, which indicates the scattered nature of the 
extension. Hence, we appear to have a monopolar or possibly 
bipolar reflection nebula emanating from S235 IRS 1 in an E-W 
direction with the blueshifted lobe on the east. Unfortunately, 
there are as yet not molecular line maps with which to look 
for the outflow. Near-IR standard imaging in molecular hydro- 
gen by M. McCaughrean (private communication) shows some 
elongated diffuse emission in the E-W direction extended for 
~ 1', which supports the picture of matter outflowing in this 
direction. 

S235IRS 4 it is associa ted with the optical nebulosity 
S235B jKrassneretai1ll982l) . The source lies ~ 1' south of 
the classical Hn region S235A. Centred on the region, there is 
also a large scale bipolar outflow oriented in the NE-SW direc- 
tion with the approaching gas towards the NE llBalrv & Ladal 
Il983t iNakano & Yoshidalll986l) . S235 IRS4 was observed in 
the H band at UKIRT. The data show no evidence of sub- 
arcsecond extended emission at a level of a 5% of the bright- 
ness peak. The source appears also unresolved in the K-band 
image from Calar Alto. Little is also known of S235IRS2, 
which appears unresolved in both images (H and K band) from 
the campaign at Calar Alto at the 5% level. 



3.9. GL2591 

This massive YSO is associated with an optical and near- 
IR nebula that extends for approximately 20" to the west 
llLenzenll 19871 lYamashita et alJI 19871 iRolph & Scarrottlll99lt 
lTamuraetalJll99lllMinchinetalJll99lh . The source illumi- 
nating the nebula (IRS1), which is inferred from the centro- 
symmetry in the optical and near-IR polarisation pattern, is 
completely obscured in the optical. The source is located near 
the centre of a large scale bipolar outflow oriented in the E- 
W direction with the blue-shifted lobe pointing to the west 
jLada et alJl984llMitchell et alJl992i) . lYamashita et ai1 ( ll987l) 
detect extended CS emission in the N-S direction that can be in- 
terpreted as a large scale torus of dense material perpendicular 
to the outflow. 

GL2591 is associa ted with a faint r adio source elong ated 
in the E-W direction (Campbell 1984a). Tamura & Yamashita 
ill 9921) also detect a group of H2 knots aligned in the E-W direc- 
tion. Most of the emission occurs towards the approaching lobe 
of the outflow, although there is also some shock-excited gas to 
the west of IRS 1 . The line of H2 knots to the E of IRS 1 appears 
to be continued by a group of H H objects aligned roughly along 
the same direction dPoetzel et alJl992l) . 

The speckle data of GL2591 taken at UKIRT show the 
presence of a faint knot 076 + 071 away at PA = 259° + 8° with 
respect to the brightness peak (figure \Y2\. The contribution to 
the flux of is 1.3% + 0.5%. This feature could be the base of the 
larger scale reflection nebulosity to the west associated with the 
blueshifted outflow lobe. 
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Fig. 11. a) Reconstructed image of S235 IRS1 in the K band taken with MAGIC during the Calar Alto campaign. The resolution 
achieved is 0727. In addition to IRS1, which is near the centre of the FOV, two other stars can be seen towards the NW and SW. 
The contours represent 0.5%, 1%, 1.5%, 2%, 3% and 5% of the peak. The noise level is ~ 0.1% although negative artifacts at a 
level of -2.5% (white in the grey-scale) appear close to the brightness (central) peak, b) The same as on panel a) but using the 
data taken during the campaign at UKIRT. In this case, the negative artifacts are at a -1% level. 
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Fig. 12. Speckle reconstructed image of GL2591 IRS1 in the 
K band at a resolution of 0721. There appears to be a faint ex- 
tension to the west, at a distance of 076 + 07 1 from the bright- 
ness peak. The contour levels are at 1.5%, 2%, 5%, 30% and 
50% of the peak. The noise level is < 1 % but the negative ar- 
tifacts (white in the grey-scale) seen to the E of the brightness 
peak are at a -1.7% level. 
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Fig. 13. Speckle reconstructed image of GL 4029 IRS 1 in the 
K band at a resolution of 0719. The contour levels are at 1%, 
5%, 10% and 50% of the peak. The noise level is < 1%. The 
grey-scale varies from -0.8% (white) to 2% (black). 
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3.10. GL4029IRS1 

This massive YSO is a member of a cluster of IR sources 
(Deharveng et al. 1997) located at a distance of 2.2 kpc 
ibecker & Fenkarfl 1 1 97 lh . GL 4029 I RS 1 is located ne ar the 
centre of a bipolar molecular outflow iSnell et alJl98«t> . Both 
outflow lobes have their peak emission coincident with the IR 
source. This maybe caused either by a low degree of colli- 
mation or by the outflow pointing towards us. GL 4029 IRS 1 
has an associated reflec ti on nebula i n the optical and near-IR 
foeharveng et alJ Il997l) . iRav et al.1 Jl990l) found several op- 
tical emission line knots emanating from GL 4029 IRS 1 that 
indicate the presence of a jet at a PA ~ 255° with a pro- 
jected ext ension of 13". Hi gh resolution observations with 
the VLA JZaoata et alJl200lT) resolve the emission at 3.6 cm 
into two components separated by 076 in the N-S direction 
with a possible counterpart to the optical jet being detected in 
the southern component. Sub-arcsecond resolution mid-IR ob- 
servations also show some extended emission emerging from 
GL 4029 IRS 1 JZavagno et aDll999t) . which appears to be as- 
sociated with the jet. 

Re cent HKL' ID speckle interferometry bv iLeinert et alJ 
J200ll) show the presence of an extended scattered light halo of 
diameter 1'.'4 + 073. The K-band speckle image taken at Calar 
Alto (figurefOl shows some signs of such an extended halo, al- 
though it is somewhat marred by artifacts from the reconstruc- 
tion procedure. Similarly, our K-band visibility function could 
be seen as being consistent with a 20% extended halo contribu- 
tion, but unfortunately the visibility function of one reference 
star with respect to the other shows a similar drop at high spa- 
tial frequencies an d so is not conclusiv e. However, it is inter- 
esting to note that ILeinert et alJ J200ll) deduced the presence 
of a patch of separate diffuse emission to the NW of the cen- 
tral source from their ID data and this is exactly where we see 
the most convincing part of the halo in our K-band image. We 
also see hints of some extended emission just to the north of 
GL 4029 IRS 1 . Neither of these correspond with the jet direc- 
tion and so their nature is unclear. 

Besides IRS1 there are two more stars, one towards the 
north at a distance of 371 ± 071 and a PA of 1.27° ± 0.03° 
and the other towards the south at a distance of 276 + 07 1 and 
a PA of 166.04 + 0.04. Two emission features are seen towards 
the NE within 1" of IRS 1, but their nature is questionable and 
neither exactly coincides withi n the positional erro rs with the 
northern radio source found bv lZapata et alJ J200ll) . 

3.11. NGC7538IRS1 

NGC 7538 IRS 1 is the brightest of a cluster of IR sources lo- 
cated in the southern edge of the optical Hn reg ion NGC 7538 
llWvnn- Williams et al.ll974[lTamura et al Jl99 lb . The cluster is 
at the centre of a bipolar molecular outflow at a PA ~ -45° that 
extends for ~ 1 pc with the blue-shifted l obe towards the NW 
jFischer et alJIl 98.4 IScoville et al JM 986t i Kameva et al J fl989t 
iDavis et al.ll998h . 

VLA observat ions at 5 and 15 GHz of NGC 7538 IRS 1 
(Campbell 1984b) show that the compact ionised gas emis- 
sion has a bipolar morphology elongated in the N-S direc- 
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Fig. 14. Speckle reconstructed image of NGC 7538 IRS 1 in the 
K band at a resolution of 0733. The contour levels are at 0.3%, 
0.5%, 1%, 2%, 25% and 50% of the peak. The noise level is ~ 
0.2%. The grey-scale varies from -0.1% (white) to 1% (black). 
In this case, the full FOV of MAGIC was used, and therefore 
the field of view is 1779 x 1779. The two parallel horizontal 
stripes that cross the image at the positions of IRS 1 and IRS2 
are artifacts from the reconstruction. 

tion (PA ~ 165 - 180°). T hey interpret this emi ssion as the re- 
sult of a collimated wind. iGaume et all {i995) found that the 
H66a recombination li nes were very broa d, indicating a wind 
speed of ~ 250 km s' MDavis et al.l (119 9 8) found shock excited 
H2 (1 - 0) S(l) emission extendin g about an arcminute north- 
wards of NGC 7538 IRS 1 (see also lBloomer et alJl998h . Bow- 
shaped features are also seen in the vicinity of IRS1, but it is 
difficult to uniquel y iden tify their powering source. 

IScoville et alJ Jl986l) found a dense core of 13 CO emission 
towards NGC 7538 IRS1 that is elongated in the E-W direction. 
The detection of a line of methanol masers that extends for 
~ 300 AU at a PA ~ 1 10° (roughly perpendicular to the outflow 
direction) may indicate the presence of an edge-on disc in this 
source (Mini eretall200ll) . 

The region shows also a large scale near-IR nebulosity 
(~ 0.3 pc) which has a monopolar shape extending towards 
the blue-shifted lobe of the la rge scale molecula r outflow with 
NGC 7538 IRS 1 at its apex jTamuraet al.lll99ll) . Its centro- 
symmetric polarisation pattern indicates that NGC 7538 IRS 1 
is the illuminating source. 

Our K band reconstructed image of NGC 7538 IRS 1 has 
a complex structure (see figure I14> . Due to the presence of 
several bright sources in the field the artifacts are particularly 
strong here, especially the two parallel horizontal stripes that 
cross the image at the positions of IRS1 and IRS2. In addi- 
tion to IRS1, IRS2 and IRS3, other near-IR sources also ap- 
pear in the 1779 x 1779 field of view of the image. The stars 
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labelle d from 'a' to 'd' also appear in the seeing-limited im- 
age of bloomer et al.l dl998t) . The stars 'e' and 'f however, 
are not seen in their image. The extended features 'A' to 'C, 
which are part of the large scale near-IR reflection nebula, ap- 
pear as regions of enh anced emission in the image presented by 
bloomer et all Jl998l) . 

It is notable that the extended emission is mainly concen- 
trated towards the N and NE of IRS1, which is roughly the 
orientation of the blue-shifted lobe in the CO outflow. Features 
'B' and "C would appear to be the result of scattered light in 
the walls a limb-brightened cavity at a PA ~ -35° with its apex 
at IRS1 and with an opening angle of ~70°. However, no bright 
sub-arcsecond nebulosity is seen that would correspond to the 
base of this large scale cavity. 

3.12. GL5180IRS1 

GL5180IRS1 is n ear the centre of a CO outflow oriented 
at a PA ~ 130° JSnell et alJ fl988). K-band imaging po- 
larimetry shows a s mall cluster wi thin extensive nebulosity 
framura etalJll99ll lYao et alJEoOol figure ^] shows our re- 
constructed image of GL5180IRS1 in the K band covering 
a field of view of 18" x 18". We find a new binary at a dis- 
tance of 1278 + 07 1 and 107.5° + 0. 1° with respect to IRS 1 . The 
two binary components (A and B in figure fT31 > are separated 
by 0/3 ±071. The component B is at a PA of 320° + 1° with 
respect the component A. Two other IR point sources are seen 
at 10/6 ± 071 and PA of 113° ± 1° (source C) and at 1278 ± 071 
and PA of 94° ±1° (source D). None of these s ources (except 
for IR S1) coincides with the IR sources found bv lTamura et alJ 
Jl99lh . No extended emission is detected. 

3.13. GL989 

GL989 (NG C 2264 IRS) is an infrared source discovered by 
lAllenl ( 1972), which is located at a distance of 750 pc. The 
presence of a molecular outflo w in the region was first de- 
scribed by Ball v & Ladal Jl983l) . The molecular outflow struc- 
ture was studied in more detail by Sch rever et alJ 1119971) . Their 
CO map appears to show that the outflow is oriented rather 
pole-on. The IR source lies at the apex of a fan-shaped neb- 
ula with an approx imate visual magnitude of 15, extending to 
the NW (Sch midtlll972|). The morphology of the nebu la was 
studied in the optical bv lScarrott & Warren-Smithl(ll989l) . They 
observed a colour gradient characteristic of reflection nebulae, 
i.e. it is redder the closer to the compact IR source. However, 
the absence of centro-symmetry around GL 989 in the polarisa- 
tion pattern, indicates that the nebulosit y is not simply a reflec- 
tion nebula illuminated by this source. Schr ever et all i 19971) 
studied the morphology of the extended nebula in the near-IR. 
In this wavelength range, the nebula is also more extended to- 
wards shorter wavelengths, which is consistent with a scattered 
origin for the emission. 

Our H-band visibility function obtained at Calar Alto in- 
dicated that some 20% of the flux was extended on scales 
larger than 1". Hints of diffuse emission to the NW are seen 
in the reconstructed image as expected if we are detecting the 
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Fig. 15. Speckle reconstructed image of GL5180IRS1 in the 
K band at a resolution of 0707. The noise level is ~ 0.5% 
although near IRS1 there are noise features with a value of 
-15%. The grey-scale varies from -1% (white) to 8% (black). 
In this case the full array was used, and therefore the field of 
view is 18" x 18". The image shows several artifacts produced 
during the reconstruction. These appear on the row through 
IRS1 and as positive and negative artifacts of the binary im- 
ages. 



larger scale nebula, but no sub-arcsecond structure was seen. 
Repeated observations with UKIRT at both H and K-bands 
failed to detect any extended emission at the 4% level. High 
resolution obser vations with NIC MOS on board the Hubble 
Space Telescope ( Thom pson et al.l 19971) also failed to detected 
any extended emission. 

3.14. NGC2024IRS2 

No nebulosity was found around this object, but a fairly bright 
companion was found 571 to the east (PA = 283°). It has a 
flux at K of about 10% that of NGC 2024 IRS2 and 20% at 
H. Hence, it is somewhat bluer than IRS2, but still a very red 
obje ct. This object is not mentioned in the mid-IR survey of 
lHaischet alJfeOOll) . 

3.15. MonR2IRS2 

This source is embedded in the Mon R2 Hn region itself and 
the strong background nebulosity caused some problems in the 
analysis. However, the source does not appear to have any sub- 
arcsecond structure. Inspection of the radially averaged visi- 
bility function yields an upper limit on any contribution to the 
total flux of about 10% level. 
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3.16. S106IR 

No sub-arcsecond structure was found in this object with up- 
per limits on a contribution to the total flux of about 3% at 
K and 6% at H. Possible reflected features are seen extend- 
ing with a wide opening angle a fe w arcseconds to the s outh 
and west in the r-, i- and z-bands bv lPersson et alJ (I1988). but 
their nature has not been confirmed polarimetrically. Recent 
near-IR Adaptive Optic s imaging of thi s source shows no 
sub-arcsecond structure (Feldtetal. 2002). Pu ga et alJ |2003) 
present a K' polarisation map of this region at a resolution 
of Of! 3. No conclusive proof of the presence of the extended 
features found by Persson et al. can be seen in this map. No 
clear polarisation pattern is observed within a few arcseconds 
ofS106IR. 



4. Discussion 

Six (GL437, MonR2IRS3, GL490, S255 IRS3, S255IRS1, 
S235 IRS1) out of the 21 massive YSOs observed show a sig- 
nificant small scale monopolar reflection nebula above the ~\% 
dynamic range of the speckle imaging presented here. In these 
cases, the nebular morphology is consistent with light being 
scattered off the dust in the walls of an outflow cavity. There is 
usually good agreement between the position angle of the sub- 
arcsecond nebula and the larger scale outflow indicators such 
as CO and larger near-IR nebulae. This agreement should be 
taken with caution though, since in some instances, large scale 
CO flows appear as multiple outflows in interferometric studies 
(e.g. iBeufher et all2003h . 

We now address the question of why only 28% of the 
sources show nebulosity at this level and the others do not, 
bearing in mind that the sample was not selected in any 
systematic way. Simulations show that by far the most im- 
portant parameter for the brightness of a reflection neb- 
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cal depth and the source is either totally invisible or only a 
faint nebula is present with no bright point source with which to 
utilise the speckle imaging technique. When the optical depth 
becomes very low, the central star will dominate the near-IR 
light such that any nebulosity will be lost in the limited dy- 
namic range of speckle imaging. Other factors such as the in- 
clination angle, the opening angle of the cavity, the density dis- 
tribution in the envelope and the dust parameters also play a 
role in determining the nebular brightness. However, these are 
secondary if the density scale is a ltered such that the line of 
sight optical depth is kept constant dAlvarez et all2004l) . 

The picture outlined above would primarily be an evo- 
lutionary one with the circumstellar matter being gradually 
cleared by the outflow, and hence there will be less dust avail- 
able for scattering. There will also be less dust obscuring the 
central star, and therefore the contrast between the reflection 
nebula and the central source will increase. At some point 
the reflection nebulae become undetectable at the 1 % level. In 
those objects where there is insufficient optical depth at K to de- 
tect the nebula, then going to shorter wavelengths should reveal 
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Fig. 16. Colour-colour diagram for the massive YSOs observed 
in our speckle campaigns. The filled circles represent the 
sources in which a sub-arcsecond conical nebula was detected. 
The open circles represent sources in which no clear extended 
emission associated with the outflow was detected. Note that 
our typical detection limit is at 1 % of the brightness peak. The 
main sequence is plotted as a grey line. The dashed lines in- 
dicat e the extin cted main sequence using the extinction law of 
iHe etalJ fl995). The arrow represents a visual extinction of 10 
magnitudes. 



the nebulae (e.g. GL490), as the scattering efficiency increases 
and the star is obscured. Unfortunately, this can often not be 
possible with speckle imaging, as the nebula becomes also too 
faint, even though it is relatively brighter than the star. 

To investigate the detection dependence on optical depth 
we have plotted our targets on a J - H versus H - K colour- 
colour diagram in figure ^] There is no clear separation of 
those with bright sub-arcsecond nebulae and those without in 
this diagram. The most we can say is that the very lightly red- 
dened sources are not showing the phenomena. This is con- 
sistent with these simply not being embedded enough to still 
produce bright nebulae relative to their central star. However, 
not all the detections are among the most reddened sources. 
Of course, using such a colour-colour diagram as an indica- 
tor of the optical depth of the circumstellar envelope is fraught 
with many complications. Some of the highly reddened non- 
detections could simply suffer from high line-of-sight extinc- 
tion that does not arise in the immediate circumstellar envelope 
of the object itself. The K-band brightness is likely to be dom- 
inated by hot circumstellar dust contributions in many cases 
whilst the shorter wavelengths can be primarily due to unre- 
solved scattered light. 

We might expect a trend with distance in that it is easier to 
detect the cavity base in the nearby objects. However, there is 
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not any indication of this in the limited sample we have studied 
here. 

Although the aim of this study was not to search for 
close companions in massive YSOs, it is worth to note that 
3 out of the 21 sources show close companions (LkHalOl, 
MonR2 IRS 3, GL961-Wa) with separations between 0717 and 
173 . These angular separations correspond to projected dis- 
tances between 150 AU and 1800 AU. This represents a bi- 
nary detection frequency of 14+8%. The statistical significance 
in our sample is fairly poor though, and more massive YSOs 
should be studied before a robust estimate on the frequency of 
binaries in massive YSOs can be achieved. In any case, it is 
worth comparing our binary detection rate with that found in 
Herb ig Ae/Be and OB sta rs. In a sample of 31 Herbig Ae/Be 
stars. iLeinert et al.l dl997l) find that ~ 16% have sub-arcsecond 
separations and ~ 29% have projected separations < 1000 AU. 
Duch ene et aD feOOlh find, in a sample of 60 OB stars, an 
18%±6% sub-arcsecond binary detection rate (projected sep- 
arations in the range 200 AU - 3000 AU). Hence, we find a 
binary detection frequency that is in good agreement with that 
found in Herbig Ae/Be and OB stars. 

5. Conclusions 

The sub-arcsecond morphology in a sample of 21 massive 
YSOs has been studied using speckle imaging in the near-IR. 
About 30% of the sources show a sub-arcsecond monopolar 
reflection nebula. These nebulae are consistent with being due 
to light scattered off the walls of the cavity evacuated by the 
blueshifted lobe of a bipolar outflow in most cases. A wide 
range of opening angles is seen and some have a gap between 
the star and the base of the nebular cone. The size of the sam- 
ple is still far from that necessary to examine trends in open- 
ing angle such as whether the cavity widen in angle in older 
sources as the outflow clears the envelope. Larger and well- 
selected samples of massive YSOs will be needed to address 
the prope rties and evolution of the outflow cavities systemati- 
cally (see iLumsden et all2002l) . 

For the majority of sources where no sub-arcsecond nebu- 
losity is seen, it could be due to them being at a slightly more 
evolved stage where their circumstellar envelope is not suffi- 
ciently optically thick. This could be checked by resolving the 
emission from the warm or cool dust in these envelopes at ei- 
ther mid-IR or sub-m illimetre wavelengths. A In a related paper 
jAlvarez et alJl2004 . we use a Monte Carlo code of the scat- 
tered light to examine in detail the detectability and the depen- 
dence of the sub-arcsecond nebulae on the physical properties 
of the circumstellar matter. 

Adaptive optics observations are capable of much higher 
dynamic ranges and will therefore be able to detect a much 
higher fraction of nebulae than was possible here with speckle 
imaging. This will also help test the scenario in which the re- 
flection nebulae fade as the outflow clears material from the 
envelope. Larger telescopes will allow the cavity to be traced 
closer to the star where collimation of the flow is likely to be 
occurring. 
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